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Relic density of neutrinos with primordial asymmetries 
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We study flavor osciUations in the early universe, assuming primordial neutrino-antineutrino asym- 
metries. Including collisions and pair processes in the kinetic equations, we not only estimate the 
degree of flavor equilibration, but for the first time also kinetic equilibration among neutrinos and 
with the ambient plasma. Typically the restrictive BBN bound on the Ve^e asymmetry indeed ap- 
plies to all flavors as claimed in the previous literature, but flne-tuned initial asymmetries always 
allow for a large surviving neutrino excess radiation that may show up in precision cosmological data. 

PACS numbers: 14.60.Pq, 98.70.Vc, 26.35.+C 



Introduction. — The state of the hot medium in the 
early universe before and around big-bang nucleosyn- 
thesis (BBN) is fixed by a small number of cosmologi- 
cal parameters, most importantly the baryon abundance 
that nowadays is measured most accurately with preci- 
sion cosmological data. In addition there are three un- 
known neutrino asymmetries, one for each flavor, that 
modify the cosmic neutrino density and, for the case of 
the I'ei^e asymmetry, influence (3 reactions of the form 
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, . primordial *He abundance [ESSE i, S, 0] 



and thus the 



It is often assumed that sphalerons equilibrate the 
lepton and baryon asymmetries, implying that all neu- 
trino asymmetries are of order 10~^ or smaller, leav- 
ing no observable trace. On the other hand, sphalerons 
have never been observed, so a cosmic neutrino asymme- 
try would raise fundamental questions about both elec- 
troweak physics and the early universe. In any event, 
improving observational information on neutrino asym- 
metries is part of the ongoing effort to determine all cos- 
mological parameters with increasing accuracy. 

Neutrino oscillations driven by the "solar" and "atmo- 
spheric" mass differences Amg^^j 



and Am^ 
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Q II are now established, whereas the third angle is 



large mixing angles sin O12 = 0.30 
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constrained by sin^ ^13 < 0.04 and could actually van- 
ish. {2a ranges, taken from Ref. In the early uni- 
verse, flavor conversions are suppressed by large mat- 
ter effects. Oscillations driven by Am^^j^^ begin only at 
T 10 MeV whereas those driven by ArTig^^j begin as 
late as 3 MeV, still early enough to achieve strong flavor 
conversions before BBN [1, [S, [U, [H. Therefore, all 
flavors must have similar asymmetries at BBN and the 
restrictive limit on the Ve degeneracy parameter, conser- 
vatively —0.04 ^ ^i,^ <; 0.07 [J], seemingly apphes to all 
flavors, unless oscillations are blocked, for instance by a 
hypothetical neutrino-majoron coupling [l^. 

Based on this reasoning it was often concluded that 



primordial neutrino asymmetries could not cause a sig- 
nificant increase of the cosmic neutrino density without 
violating the BBN bound on . Assuming that neu- 
trinos indeed reach perfect kinetic and chemical equilib- 
rium, the possible increase of the cosmic radiation den- 
sity is small, AiVeff = f + f {^./nf < 0.006. 
Here, the usual "effective number of neutrino families" 
measures the cosmic radiation density, after e+e^ anni- 
hilation, by prad/p^ = 1 + (7/8) (4/11)^/3 iVes . 

We stress, however, that Refs. , [ill, [l2| do not claim 
that neutrinos after flavor conversions follow Fermi-Dirac 
distributions in terms of a common degeneracy parameter 
and the temperature of the ambient plasma. Refs. 11. 1^ 
showed that, if Oi^ — and if one ignores collisions, syn- 
chronized MSW-like flavor conversions driven by Am'^^^ 
and 9i2 achieve comparable vi? asymmetries among all 
flavors, but this is not the same as kinetic and chemical 
equilibrium. Ref. [lo| included collisions in a schematic 
way, showing that the coherence among flavors was ef- 
ficiently destroyed. However, the transfer of entropy to 
the electromagnetic plasma was not considered and the 
final neutrino distributions not studied. 

If all initial asymmetries have the same sign, oscil- 
lations alone are enough to achieve approximate flavor 
equipartition and the BBN limits on the Vgi^e asymmetry 
will apply to all flavors. A dangerous case arises, how- 
ever, if one begins with opposite asymmetries in two fla- 
vors. Oscillations, with or without collisions, will achieve 
approximate flavor equipartition separately among neu- 
trinos and antineutrinos, so in the end all asymmetries 
vanish. On the other hand, the original excess radiation 
remains. The efficiency of viy — s- e^e" entropy transfer 
determines if the BBN limits on the Vei^e asymmetry still 
preclude a large ATVeff. The main point of our paper 
is to study this case that was not covered in previous 
treatments [l3,[ni[il- 

Numerical Approach. — We describe neutrino distribu- 
tions by 3 X 3 matrices in flavor space Qp for each mode. 
The diagonal elements are the usual occupation numbers 
whereas the off-diagonal ones encode phase information. 
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We need to solve the equations of motion (EOMs) [lj,lla| asymmetry, is 
^Qp^[^p,Qp]+C[gp,gp] (1) 



and similar for the antineutrino matrices Pp. The first 
term on the r.h.s. describes oscillations with 



12C(3) 



(3) 



2p \ 3m;;, 



(2) 



where p — |p|, M is the neutrino mass matrix, g — 
J gp d'^p/(27r)'^ and similar for g. The small baryon den- 
sity implies that the dominant matter term in the os- 
cillation matrix is equal for neutrinos and antineutrinos, 
where E is the 3x3 flavor matrix of charged-lepton en- 
ergy densities [3, [l6| . The matrix D.p for antineutrinos 
is the same with M^/2p — M^/2p. The term propor- 
tional to p — p is responsible for synchronized oscilla- 
tions 

Solving these EOMs with the full collision terms is nu- 
merically challenging and not necessary for a first ex- 
ploratory study. Collisions have two main effects: They 
break the phase coherence of mixed flavor states and they 
create or destroy neutrino pairs. The first effect is effi- 
ciently implemented with p-dependent damping factors 
for those gp elements that are off-diagonal in the weak- 
interaction basis. Here we follow exactly the earlier work 
of Dolgov et al. [lo| in that we use the damping factors 
spelled out in Eqs. (60) and (61) of Ref. Essentially 
one assumes that neutrinos collapse into weak-interaction 
eigenstates each time they interact in a process that dis- 
tinguishes among flavors [17 1. 

The crucial new ingredients are collisions and pair pro- 
cesses for the diagonal gp elements, allowing the neu- 
trino distributions to achieve equilibrium with the ambi- 
ent plasma. With small modifications we here use an ex- 
isting numerical code JSj where the collision integrals are 
solved without approximation, using the same approach 
as Ref. [l^. For example, we obtain exact Fermi-Dirac 
distributions if the relaxation time is large enough. 

So in our implementation the diagonal elements change 
by collisions and oscillations, whereas the off-diagonal el- 
ements change by oscillations and damping. We have 
checked that the damping factors of Ref. [15[ and the col- 
lision rates used in the code [13] are consistent with each 
other. Changing the damping factors by as much as a fac- 
tor of 2 has only a minimal impact on the final AiVeff . Of 
course, our combined treatment of oscillations, damping 
of flavor coherence, and collisions and pair processes is 
only approximate and would be questionable if one were 
to aim at a precision determination of AN^ff. However, 
our relatively simple extension of existing tools provides 
a first reasonable estimate of the surviving AA^eff- 

Opposite Equal Asymmetries. — Our main example is 
the "dangerous case" of two large but opposite initial 
asymmetries. To be specific, we first consider an exactly 
vanishing global lepton asymmetry. The initial asymme- 
try in one flavor z/q, defined in analogy to the baryon 



A vanishing global asymmetry requires for the initial de- 
generacy parameters 2 {n'^^^^ + ^tj + (•^^C^e + = 0. 

In Fig. [1] we show the evolution of p^/p^, properly 
normalized so that it corresponds to A^cff at early and 
late times, first for the baseline case without asymmetries 
(dotted hues). The fast drop of p^/p^ a.t T ^ 0.2 MeV 
represents photon heating by e+e^ annihilation. At late 
times the dotted lines end at Nes = 3.046 instead of 3 
because of residual neutrino heating [l3| • 

Next we show (solid lines) the evolution for our main 
example where = —0.64 and = 1.17, correspond- 
ing to an initial value iVcff — 4.0. The neutrino mass dif- 
ferences and 023 are fixed to their best-fit values. In the 
upper panel we assume Oi^ — and for the upper curve 
also 6*12 = (no oscillations), leading to an unchanged 
final A^off = 4.0. The bottom curve is for sin^ 612 = 0.5 
(maximal mixing) , showing that at T-y ~ 3 MeV the neu- 
trino density depletes relative to photons, when oscilla- 
tions driven by Ato^^j begin and the fiavour asymmetries 
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FIG. 1: Evolution of the neutrino energy density. The vertical 
axis is marked with A'^off , left before e^e~ annihilation, right 
afterwards. Dotted lines: Baseline case without asymmetries. 
Top panel: 6ii = 0, the solid lines correspond to sin^ 612 — 0, 
0.1, 0.3, and 0.5 from top to bottom. Bottom panel: sin^ 612 = 
0.3, the top solid curve is sin^ 613 — 0, while the other solid 
curves correspond to sin^ 613 — 0.04 with normal (top) and 
inverted (bottom) hierarchy. 
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p/T 

FIG. 2: Final and Ve spectra for ^13 — and sin^ S12 = 0.3. 
The dotted lines are Fermi-Dirac spectra producing the same 
VcVe asymmetry and the same energy density pi,^ + pa^ . 



are reduced. The main effect of collisions is to keep the 
neutrino spectra closer to their equilibrium form: Phys- 
ically the excess of entropy is transferred from neutri- 
nos to the electromagnetic plasma, cooling the former 
and heating the latter. Practically the same behavior 
is found for the next curve representing sin^ 612 = 0.3, 
the experimental best fit, where in the end TVeff = 3.24. 
A further curve is shown for the value sin^ ^12 — 0.1, 
remaining far from equilibrium until the end. In the bot- 
tom panel we use the measured sin^ 9i2 = 0.3 and repeat 
the case ^13 = (top curve). The other two curves are 
for sin^ Oxz = 0.04, near the experimental limit, for the 
normal and inverted hierarchies. 

As expected, a nonvanishing ^13 together with Am^tj„ 
drives the system towards equilibrium at a time when 
collisions are still efficient. In our numerical example, a 
significant deviation from equilibrium survives only for 
a very small ^13. In this case the outcome depends on 
6*12 as found in the original studies 10, HI, [l2|- A small 
9i2 does not lead anywhere close to flavor equilibrium, a 
large value leads to near-equilibrium. 

For 6*13 = and sin^ 6*12 — 0.3 we show the final v^, 
and Df. energy spectra in Fig. O We compare them with 
Fermi-Dirac spectra that produce the same Vgi^e asym- 
metry and the same energy density p^^ -\- p^^ , implying 
f^y^ — 0.185 and a T^^ only 5% higher relative to the equi- 
librium value. These spectra are much closer to kinetic 
equilibrium than expected. 

Asymmetry Scan. — Next we study the final deviation 
from equilibrium for a range of initial asymmetries < 
Ci/c ^ 3, assuming a vanishing global asymmetry as be- 
fore. In Fig. [3] we show the final outcome in terms of the 
residual excess radiation density AiVofi and and T^^ 
chosen for a Fermi-Dirac distribution that approximates 
the final v^^ and Ve spectra such that we reproduce the 
numerical asymmetry and energy density. We show these 
results for the experimental best-fit value sin^ 612 = 0.3 




FIG. 3: Parameters for the final Ue and i/e spectra as well as 
the final AA'off as a function of the initial . The solid lines 
correspond to 6113 = 0, the dotted lines to sin^ 6113 = 0.04. 
In the bottom panel, the dashed line is the surviving AA^cff 
when the final is in the range allowed by BBN. 

and for 6*13 = as well as sin^ ^13 — 0.04. 

Qualitatively we confirm that kinetic equilibrium is 
rather good in that the final neutrino temperature is close 
to what it should be in equilibrium. Moreover, we also 
find that the value of ^13 is crucial in determining how 
close the system comes to equilibrium, which however is 
never perfect for any value of 6*13. If sin^ < 10""^, the 
outcome is very close to a vanishing ^13. 

We note that for large initial the final electron neu- 
trino degeneracy can be so large that it is excluded by 
BBN. However, we are not bound to the assumption of 
an initially vanishing global asymmetry, but we can use 
slightly more negative initial ^i^^ values such that the fi- 
nal is pulled into the range allowed by BBN, yet a 
significant excess radiation density survives as shown in 
in Fig. [21 being quite similar to ANes for the simple case 
of a vanishing global asymmetry. In other words, up to 
small adjustments of initial conditions, the simple case 
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of a vanishing global asymmetry gives us a fair idea of a 
possible surviving residual radiation density. 

No Initial v^Ve Asymmetry. — Another initial configu- 
ration for a vanishing global asymmetry is = and 
^i/^ = —^1/^. If the latter are exactly equal but op- 
posite, v^tVT oscillations driven by Am^^jj^ and 023 are 
suppressed by neutrino-neutrino interactions in that the 
synchronized oscillation frequency vanishes exactly [lo| . 
However, oscillations driven by Am^Qi and Q\i achieve 
approximate equilibrium, the transformations beginning 
at T '--^ 3 MeV as before. We have not studied this case 
in detail because it is numerically more challenging. We 
note that "blocking" the v^Vr oscillations requires a very 
precise cancellation and thus a more special case than 
before where an exactly vanishing global asymmetry was 
only a choice of convenience. 

Conclusions. — The first studies of primordial neutrino 
flavor equilibration were concerned with the role of 
the "solar" neutrino mixing parameters and found that 
approximate global flavor equilibrium is only achieved 
for the now-established "large mixing angle solution" 
[lol . 11, [l2|. While this conclusion is correct, we here 
stress that it must not be over interpreted. Oscilla- 
tions driven by the solar parameters take place so late 
(T ~ 3 MeV) that neither perfect flavor equilibrium nor 
thermal equilibrium with the ambient plasma is achieved, 
although a large value for 9i3 helps to come closer because 
flavor transformations begin earlier. (If sin^ ^13 ^ 10"'^ 
it plays no signiflcant role.) For 6*13 = we have con- 
structed an explicit example with nearly equal but oppo- 
site asymmetries between ^'e and the other flavors such 
that the BBN limit on is satisfied, yet a large ATVoff 
survives. 

Precision cosmology has advanced in such strides that 
today one can derive nontrivial bounds on the cosmic 
radiation density that however still allow for AiVcff of a 
few 0, El, [13, Sin, 111. It is conceivable that in future 
one can do much better and constrain or detect a value of 
AiVeff significantly smaller than unity [1^, [l^, [H, [H, 113] , 
for instance with very precise CMB data from Planck. 
Such results remain complementary to BBN limits on 
^1,^ . If precision cosmology were to turn up a nonva- 
nishing AN^s, it could be a remnant of a primordial 
neutrino asymmetry, although the exact interpretation 
would strongly depend on 613. In this sense, future mea- 
surements of 013 are important for the interpretation of 
cosmological precision data. 

In summary, the established neutrino mixing parame- 
ters together with BBN constraints on typically as- 
sure that the surviving excess radiation density is im- 
measurably small, in agreement with the previous liter- 
ature. However, imperfect kinetic and chemical equilib- 
rium achieved by the solar mixing parameters allows for 
exceptions. In particular, large but opposite primordial 
asymmetries of two fiavors can be fine-tuned such that 
the limit is respected, yet a large excess radiation 



density survives with AiVogf of order unity or larger that 
may be detectable in future cosmological precision data. 
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